A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
13.25.2 11.19.1 10.09.1 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



BeppoSAX unveils the nuclear component in NGC 6240 

p. Vignatii, S. Molendi\ G. Matt^, M. Guainazzi^, L.A. Antonelli'', L. Bassani^ W.N. Brandt^ A.C. 
Fabian^, K. Iwasawa^, R. Maiolino®, G. Malaguti^, A. Marconi*, and G.C. Perola^ 

^ Istituto di Fisica Cosmica "G.Occhialini" , Via Bassini 15, 1-20133 Milano, Italy 

Dipartimento di Fisica "E. Amaldi", Universita degli Studi Roma Tre, Via della Vasca Navale 84, 1-00146 Roma, Italy 
^ Astrophysics Division, Space Science Department of ESA, ESTEC, Postbus 299, NL-2200 AG Noordwijk, The Netherlands 

* Osservatorio Astronomico di Roma, Via dell'Osservatorio, 1-00044 Monteporzio-Catone, Italy 

^ Istituto Tecnologie e Studio Radiazioni Extraterrestri, CNR, Via Gobetti 101, 1-40129 Bologna, Italy 

^ Dept. of Astronomy and Astrophysics, The Pennsylvania State University, 525 Davey Lab, University Park, PA 16802, U.S.A. 

^ Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, U.K. 

* Osservatorio Astrofisico di Arcetri, L. E.Fermi 5, 1-50125 Firenze, Italy 

To appear in A&A Letters 



' Abstract. The direct nuclear X-ray emission of the ultra- 

. luminous IRAS galaxy NGC 6240 has been observed for 

' the first time by BeppoSAX.. It is seen through an absorber 

! with A^H 2 X 10^4 cm-2 rj.^^ 2-10 keV X-ray luminosity 

' of > lO"*^* erg s~^ definitely proves that a powerful AGN is 

. hosted in this galaxy, despite the LINER optical appear- 
ance. Assuming the typical Lx/^boi ratio for quasars, our 

. result implies that it is the AGN and not the starburst 

' that dominates the energy output of NGC 6240. 

Key words: X-rays: galaxies - Galaxies: Seyfert - Galax- 
ies: individual: NGC 6240 



1. Introduction 

One of the major open problems with Ultraluminous In- 
frared Galaxies (ULIRG) is whether they host active nu- 
clei at their centers. In at least a fraction of objects, there 
is circumstantial evidence that this is the case. One of the 
best examples in this respect is NGC 624Cu, a well known 
source in which the presence of two nuclei (Fried & Schulz 
1983) suggests an ongoing galaxy merger. Indications of 
the presence of an AGN come mainly from X-rays: the 
flat 2-10 keV X-ray spectrum and the prominent iron line 
complex observed by ASCA (Iwasawa & Comastri 1998) 
are very similar to those of NGC 1068 (Iwasawa, Fabian 
& Matt 1997; Matt et al. 1997 and references therein) and 
suggest reflection from cold and warm matter of an oth- 
erwise invisible nuclear component. On the other hand, in 

^ If LiR > 10^^ Lq is the limiting criterion for ULIRG, strictly 
speaking NGC 6240 belongs to this class only if Ho <65 km 
s"^ Mpc"^ Genzel et al. (1998). 



the optical spectrum all the diagnostic line ratios point 
towards a classification as a LINER (e.g. Veilleux et al. 
1995), and the ISO SWS diagnostic diagram places NGC 
6240 in the region of star formation dominance (Genzel et 
al. 1998). In soft X-rays the spectrum is fairly complex, 
and best modeled by multi-temperature plasma emission, 
very likely associated, at least partly, with the powerful 
starburst and related wind present in this galaxy, a hy- 
pothesis strengthened by the ROSAT HRI discovery of 
rather luminous (^10"*^ erg cm^^ s^^) emission extended 
over tens of kpc (Komossa, Schulz & Greiner 1998; Schulz 
et al. 1998). On the other hand, the point-like soft X- 
ray component is possibly due to scattering of the nuclear 
radiation (Komossa et al. 1998). 

A direct measurement of the nuclear luminosity is 
needed to assess the relative importance of AGN and star- 
burst emission. The only way to do that is to explore 
the hard X-ray band, which may reveal the direct nu- 
clear component. If the absorbing matter is moderately 
thick to Compton scattering (i.e. iVn ~a fewxlO^'' cm~^), 
the photons emerge above 10 keV, where the most sensi- 
tive instrument presently available is the PDS onboard 
BeppoSAX. This has already added two sources (Circinus 
Galaxy, Matt et al. 1999; Mrk 3, Cappi et al. 1999) to the 
small list of "moderately-thick" Seyfert 2s, which was pre- 
viously composed of only one object, NGC 4945 (Iwasawa 
et al. 1993; Done et al. 1996). 

2. Observation and data reduction 

BeppoSAX (Boella et al. 1997) observed NGC 6240 on 
1998 August 14-17. In this paper, data from three instru- 
ments are discussed: the LEGS (0.1-10 keV), the MEGS 
(1.5-10 keV), and the PDS (13-300 keV). 
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Data reduction and analysis follow the standard cri- 
teria, as described in Matt et al. (1997). The total ef- 
fective exposure time was 52.5 ks in the LEGS, 119.4 ks 
in the MEGS and 114.4 ks in the PDS. We accumulated 
and examined images for the MEGS instruments. Gom- 
paring the radial profiles in three energy ranges (1.5-4, 
4-7.5, 7.5-10 keV) with the corresponding Point Spread 
Function, no evidence for extended emission is found. 
Light curves and spectra were extracted within regions 
of 4' and of 2' radius centered on the source for MEGS 
and LEGS, respectively (the small LEGS extraction ra- 
dius is chosen to avoid contamination from a nearby, soft 
source, not visible in the MEGS). The light curves are 
consistent with being constant. The resulting count rate is 
(8.46±0.67)xl0-3 counts s'^ (0.1-4.5 keV) for the LEGS 
and (2.46±0.07)xl0-2 counts s^^ (1.5-10 keV) in the 
MEGS. No significant variability was found in the PDS 
light-curve either, with a count rate of 0.38±0.04 counts 
s^^. We restricted the spectral analysis to the 0.1-4.5 
keV and 1.5-10 keV energy bands for LEGS and MEGS 
respectively, where the latest released (September 1997) 
response matrices are best calibrated. Standard blank- 
sky files provided by the _BeppoSAX Science Data Genter 
(SDG) were used for the background subtraction. 

3. Spectral analysis and results 

In Fig. 1 we compare the >2 keV spectra of NGG 6240 
and NGG 1068. Both spectra have been normalized to the 
Grab spectrum; since the latter is a power law with pho- 
ton index ~ 2, we are basically showing the i/F^ spec- 
trum. The MEGS spectra of the two sources are similar 
in shape, not surprisingly as we know that in this energy 
range both sources have been successfully fitted by two 
reflection components, warm and cold (e.g. Iwasawa et al. 
1997; Iwasawa & Gomastri 1998). While below 10 keV 
the flux of NGG 6240 is the lowest, the opposite occurs 
in the PDS band. Glearly, in NGG 6240 a new component 
emerges at high energies. No known bright sources are 
present in the field of view of the PDS, and the probabil- 
ity of a serendipitous source in the 1.3° PDS field with a 
flux equal or larger than that of NGG 6240 is < 10"^, if 
the 2-10 keV ASGA LogN-LogS (Gagnoni et al. 1998) is 
adopted, and a power law spectrum with photon index 1.8 
is assumed. Thus the best explanation for this component 
is nuclear emission piercing through an absorber with col- 
umn density Nh a fewxlO^^ cm i.e. Gompton-thick 
but still permitting partial transmission above ^10 keV. 
In Fig. 2, the broad band spectrum is fitted following Iwa- 
sawa & Gomastri (1998): an excess is evident in the PDS 
band. We therefore modeled the continuum above 2 keV 
with: i) a power law, representing the intrinsic emission, 
with a photon index fixed at a typical value of r=1.8, 
absorbed by a column density allowed to be free (leav- 
ing the index as a free parameter, a value consistent with 
the adopted one, but very loosely constrained, is found); 



ii) a second, unabsorbed power law with the same index, 
standing for the warm reflection component; in) a pure 
cold reflection component (PEXRAV model in XSPEG) 
generated by a power law with the same index. 
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Fig. 1. Comparison of the MECS and PDS spectra of NGC 
1068 (open stars) with those of NGC 6240 (dots). Both spec- 
tra have been normalized to the Crab spectrum, which is a 
power-law with a photon index of ~ 2. 

As we are interested here in the high energies, and be- 
cause below 2 keV ASGA is superior to BeppoSAX as far 
as both sensitivity and spectral resolution are concerned, 
we content ourselves with adopting in this band the best 
fitting model of Iwasawa & Gomastri (1998). It consists 
of a multi-temperature optically thin plasma, which prob- 
ably originates in a powerful starburst. In particular we 
use a two-temperature MEKAL model, where the absorp- 
tion of the coldest MEKAL was fixed at the Galactic value 
(Nh = 5.8 X 10^° cm^'^), while the absorption of the other 
MEKAL was allowed to be a free parameter. We applied 
the latter absorption to all the AGN components. The 
envisaged scenario is then one in which an external star- 
burst region, associated with the cooler MEKAL compo- 
nent, absorbs emission from the inner regions of NGG 6240 
(e.g. Fabian el al. 1998). 

A prominent and broad iron K line complex, very sim- 
ilar to that seen in NGG 1068, is also evident. Because of 
this similarity we assume that the line is actually a blend 
of: the 6.4 keV Ka and 7.06 keV K/3 fluorescent lines from 
cold iron (the intensity of the latter has been forced to be 
equal to 1/9 that of the former, as expected from atomic 
physics); the 6.7 and 6.97 Ka recombination lines from 
He- and H-like iron. In summary the complete model, F, 
we used to fit our data can be written as: 

F = Ag{Asb[At{PL) + Rw + Rg + Mh + 4GL] + Afc}, 

where: Aq is the absorption associated with the Galactic 
column, AgB is the absorption related to the starburst. 
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At is the absorption acting on the nuclear emission (in- 
cluding the Compton cross section), possibly associated 
with the molecular torus shrouding the nucleus, PL is the 
power-law modeling the nuclear component, i?w and -Rc 
are respectively the warm (optically thin) and cold (op- 
tically thick) reflection components, Mh and Mc are re- 
spectively the cold and hot thermal components and AGL 
are the four Gaussian lines modeling the iron blend. 



Table 1. Best fit parameters, see text for details. 
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towards other directions are now scattered into the line of 
sight, and the nuclear luminosity required to explain the 
observed flux is therefore lower. 

The same fit but without Rc is also acceptable 
(X^ = 133/120 dof). In this case, the observed 6.4 keV 
iron line would entirely originate in the thick absorbing 
medium, which is possible provided that the covering fac- 
tor is not much smaller than unity, as we have verified 
by Monte Carlo simulations. Therefore, the presence of 
the cold reflection component is not certain in this source. 
If i?w is instead excluded, the fit is significantly worse 
(X^ = 139/120 dof). 
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Fig. 2. LEGS, MEGS and PDS spectra of NGC 6240 and 
best fitting model and residuals, when the Iwasawa & Gomas- 
tri (1998) model is adopted. A further component is clearly 
emerging in the PDS band. 



In Fig. 3 we show the deconvolved best fitting model 

with all the components. The best fit parameters arc listed 
in Table 1. The fit is acceptable, xVd-o.f.=131/119. If 
we do not include the transmission component, and try 
therefore to account for the spectral component emerg- 
ing in the PDS band with pure Compton reflection (as in 
NGC1068, Matt et al. 1997), the fit is significantly worse 
(X^ = 179/121 dof; sec Fig. 2). It must be noted that the 
absorber of the nuclear emission. At, is Compton-thick. 
A proper model, including Compton scattering within the 
matter (see Matt et al. 1999 and Matt, Pompilio & La 
Franca 1999 for details on this model) should have then 
been used, if this absorber has a significant covering factor 
(as for example for a torus geometry with no significant 
vertical density gradient). When this is done, however, 
the parameters of the various components do not differ 
significantly from those obtained with the simpler model, 
apart from the flux associated with the intrinsic emission, 
which turns out to be about one third that reported in 
Table 1. This is because a fraction of the photons emitted 



4. Discussion 

The main result of the BeppoSAX observation of 

NGC 6240 is the discovery and measurement of the nu- 
clear emission. The 2-10 keV nuclear luminosity, after 
correction for absorption, is 3.6x10^^ erg (assuming 
Hq=50 km s~^ Mpc^^), or about one third of this if the 
absorber covers a large solid angle, like in a torus. 

As discussed above, the presence of cold reflection is 
not unambiguously established in this source. Let us, how- 
ever, assume that it does exist, and take the best fitting 
parameters listed in Table 1. The ratio between the 2-10 
keV cold reflection and intrinsic fluxes is about 0.3% for 
the higher value of the intrinsic luminosity and 3 times 
larger for the lower value. Assuming a torus geometry like 
that adopted by Ghisellini, Haardt & Matt (1994), and 
consistently the lower flgure for the intrinsic luminosity, 
we derive an inclination angle of the torus axis with re- 
spect to the observer of ~35° (~45° if the upper figure 
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LINER that, when observed in X-rays, turns out to be an 
AGN (see e.g., NGC 4945, Done et al. 1996; NGC 1052, 
Guainazzi & AntoneUi 1999). Differently from them (but 
similarly to two other ULIRG, IRAS 23060+0505: Brandt 
et al. 1997 and IRAS 20460+1925: Ogasaka et al. 1997), 
NGC 6240 has a high X-ray luminosity, and it may be 
called an obscured QSO. While it would be premature to 
draw conclusions on ULIRGs as a class, one cannot help 
noting that the so far elusive type 2 QSO should be better 
searched for in X-rays rather than in the optical band. 
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Fig. 3. Deconvolved best fitting model for NGC 6240. The 
thick line represents the total model, the thin lines the spectral 
components, labeled as in Table 1. 



is instead adopted). It is worth noting that both NGC 
6240 and Circinus Galaxy (Guainazzi et al. 1999) appear 
to be observed at moderate inclinations, while NGC 1068, 
which is totally obscured, is possibly seen almost edge-on 
(Greenhill et al. 1996, Matt et al. 1997; but see Kishimoto 
1999 for a different point of view). This would suggest a 
decrease in the density of the X-ray absorbing medium 
with increasing distance from the equatorial plane. 

The ratio between 2-10 keV warm reflection and in- 
trinsic fluxes is about 0.5% (1.5% for the lower luminos- 
ity). The He- and H-like lines have EW with respect to 
the warm reflection component of about 1 and 0.8 keV, 
respectively. From these values we estimate the column 
density of the warm mirror (e.g. Matt, Fabian & Brandt 
1996; Guainazzi et al. 1999), A'^H.warm, to be between 10^^ 
and 10^^ atoms cm~^, roughly consistent with the value 
estimated from the ratio between the normalizations of 
the warm reflection and of the intrinsic emission. 

The 2-10 keV X-ray luminosity is 10 to 30 times 
lower than the IR luminosity (calculated as prescribed by 
Mulchaey et al. 1994), consistent with the usual values for 
Seyferts (e.g. Mulchaey et al. 1994; Awaki 1997). In the 
case of NGC 6240, therefore, a ULIRG turns out by no 
means to be a peculiar AGN, as far as X~ray properties 
are concerned. If L2-10 kcv/^boi ~0.03, as typical for QSO 
(Elvis et al. 1994), our result implies that it is the AGN 
that dominates the energy output, and not the starburst 
as deduced by Genzel et al. (1998) on the basis of the ISO 
spectrum. 

In the optical, NGC 6240 is classified as a LINER. 
Possibly, the NLR is obscured by the starburst region, 
whose column density, from our analysis, corresponds to 
Av -^10. It is worth noting that NGC 6240 is not the first 
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